Thermopower induced by a supercurrent in superconductor-normal-metal structures 
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We examine the thermopower Q of a mesoscopic normal-metal (N) wire in contact to supercon- 
ducting (S) segments and show that even with electron-hole symmetry, Q may become finite due 
to the presence of supercurrents. Moreover, we show how the dominant part of Q can be directly 
related to the equilibrium supercurrents in the structure. In general, a finite thermopower appears 
both between the N reservoirs and the superconductors, and between the N reservoirs themselves. 
The latter, however, strongly depends on the geometrical symmetry of the structure. 

PACS numbers: 74.25.Fy, 73.23.-b, 74.45.+C 



Thermoelectric effects in electrical conductors typi- 
cally result from the asymmetry of the Fermi sea between 
the electron (E > Ep) and hole-like (E < Ep) quasipar- 
ticles. This is illustrated by the Mott relation |l| for the 
thermopower 
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This relates the potential difference AV generated by 
the temperature difference AT to the energy dependence 
of the conductivity a due to the asymmetry above and 
below the Fermi sea. 

In metals, the electron-hole asymmetry is governed by 
the parameter ksT / ' Ef arising from the next-to-leading 
term in the Sommerfeld expansion. At sub-Kelvin tem- 
peratures, this leads to a very small Q, typically below 
10 nV/K. However, recent experiments |2j, Li-Li measur- 
ing the thermopower in normal-metal wires connected to 
superconducting electrodes indicate that it exceeds this 
prediction at least by an order of magnitude, and, more- 
over, show that Q oscillates with the phase difference 
between the two superconducting contacts. 

Mott relation is expected to fail in the presence of 
the superconducting proximity effect when the geomet- 
rical symmetry in the measured sample is broken ||. 
Our aim is to show that with nearby superconductors, 
normal-metal circuits can show a thermoelectric effect 
independent of electron-hole-symmetries, since the prox- 
imity effect couples the temperatures to the potentials 
through the supercurrent. This effect is at least two or- 
ders of magnitude larger than that predicted by Eq. (JJJ 
(c.f. Fig. D}. 

We discuss the system shown in the inset of Fig. [2 with 
a supercurrent flowing between the two superconducting 
elements. Our main result (Eq. JSJ) states that the ther- 
mopower Qns between the N and S parts of the structure 
is proportional to the difference in the supercurrents at 
the temperatures T%, T 2 of the two N electrodes. More- 
over, we obtain a similar result (Eq. iJTDJ) for the ther- 
mopower Qnn between the normal parts in a geometri- 
cally asymmetric structure. This can be understood phc- 
nomenologically as follows. If T\ / T2, the temperature- 
dependent equilibrium supercurrent Is(T\) in wire 3 



is different from Is{T2) in wire 4. (For this qualitative 
picture, we approximate these wires to be at the temper- 
atures Ti, T2O Thus, a compensating effect must arise 
to guarantee the conservation of currents. Should the N 
reservoirs be kept at the same potential as the supercon- 
ductors, a quasiparticle current / qp oc (Js(Ti) — Is(T 2 )) 
from them to the superconductors would balance the dif- 
ference. However, when no current is allowed to flow in 
wires 1 and 2, a compensating N-S potential difference 
Vn-Vsou R{Is{Ti) - Is(T 2 )) is induced instead. The 
induced potentials oscillate with phase differences, sim- 
ilarly to the supercurrent, and may differ in the two N 
reservoirs, especially in asymmetric structures. 

In the following, we concentrate on the diffusive limit, 
and model the setup with the Keldysh-Usadel equa- 
tions Q. These equations assume electron- hole symme- 
try, under which Eq. JJJ predicts a vanishing Q. 

The Keldysh-Usadel equations are formulated in terms 
of the quasiclassical Green's functions (which are matri- 
ces in the Keldysh-Nambu space), but here we use their 
0-parametrization Q, for convenience. It reduces the 
problem to two sets of equations, the spectral equations 
and the kinetic equations. In terms of the parameteriz- 
ing functions x and the phase and the proximity effect 
strength, the spectral equations in normal metals are 

DV 2 = -2iEsmh{8) + \d{Vx) 2 sinh(26>) , (2a) 
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Here, the factor D is the diffusion constant of the normal 
metal, and E is the energy with respect to the super- 
conductor potential. The kinetic equations are expressed 
using the symmetric and antisymmetric parts, Jt and 



fh lij- of the electron distribution function f(E, 
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The equations imply that the spectral current densities 
jl and jr are conserved (we neglect inelastic scattering). 
Thus the observable charge and energy current densities, 



3c 



2e 



jrdE and jq = 



2e 2 



Ej L dE, (4) 



2 



are also conserved. 

The coefficients T>l, T^t-, ^"and js appearing in the ki- 
netic equations are obtained from the spectral equations: 



V L/T = - (1 + |cosh6»| 2 =f |sinh6»| 2 cosh(2Im[x])) , (5a) 

T= i| sinh#| 2 sinh(2Im[ X ]), j s = Im\j E ]. (5b) 

Here, T>l and T>x are the local spectral energy and charge 
conductivities, and js is the spectral density of the super- 
current 0. The factor T arises in the formalism and has 
an effect on the thermopower. The normal-state values 
of these coefficients are T>l = T>x = 1 and js = T= 0. 

At nodes of wires, assuming clean metallic contacts, 
the functions 9, x an d / are continuous and Kirchoff- 
like "spectral current conservation laws" 9] imply that 
Aa N j L / T and Acr N js are conserved. Here, A is the 
cross-sectional area of a wire and <j N the normal-state 
conductivity. At clean metallic reservoir contacts most 
of the functions get their bulk values 0. However, for 
energies below the superconducting energy gap A, the 
valid boundary conditions at superconductor interfaces 
are j'l = 0, prohibiting the energy flow, and fx = 0, 
assuming no charge imbalance in the superconductors. 

The coefficients js and Tcouple the energy and charge 
currents together, and give rise to a finite thermopower. 
Moreover, these coefficients oscillate with the phase dif- 
ference in the system, and thus the value of the ther- 
mopower should also oscillate. When there is no phase 
difference, T— js — 0, and the thermopower vanishes. 

The energy scale of temperatures and potentials is 
specified by the Thouless energy 
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corresponding to a wire of length L y|. Moreover, Ex of 
the link between the superconductors is a natural energy 
scale for the spectral equations. As long as Ex <C A, the 
results can be scaled to fit all systems with similar ratios 
of wire lengths and areas. 

Since there are no general analytical solutions to the 
problem, we solve the spectral equations numerically, and 
make a few approximations to solve the kinetic equations. 
However, the data shown in the figures is obtained nu- 
merically without any approximations. 

First, we note that the "local potential" fx is generally 
small (as shown by the numerical results), as are the in- 
duced potentials at the reservoirs. Thus, we can neglect 
the terms proportional to it in the kinetic equation (|3aH . 
Physically this means that we mainly neglect the effect 
of supercurrent on the energy currents and the temper- 
atures. (If the potentials were large, the omitted term 
would be the source for a Peltier-like effect [13.) With 
this approximation, we integrate the kinetic equations, 



which yields the connection between the spectral current 
densities and the distribution functions fh(x) and fx(x) 
at the ends of a wire of length L: 

jL = j^(f L (L)-f L (0)), (7a) 
Jt = (ML) - fx(0))+3sf L (0) (7b) 
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Here M L / T = \ f Q L drr are the dimensionless spec- 
tral energy and charge resistances. To simplify the final 
result, we also approximate Pi = 1 in 01, and T>x = 1 
in the latter term in (|7b|) . since the variation in T>l (away 
from superconductor interfaces) and T>x with respect to 
the energy is smaller than that of the other coefficients. 
Numerical results verify that this does not affect the re- 
sult crucially. The energy-dependent 1/Mx as a coeffi- 
cient for fx causes an important temperature dependence 
of the conductance, so we retain it. 

Using Eqs. J7J and the conservation of spectral cur- 
rents, we obtain a linear system of equations for the spec- 
tral current densities. They can be solved with respect to 
the given temperatures and potentials in the reservoirs, 
with different results for \E\ < A and \E\ > A, due to the 
different boundary conditions. Next, we integrate over 
the energy to obtain the observable current densities, af- 
ter which we require the condition j Cj \ — j c ,2 — 0. To 
solve the resulting equations for the small induced poten- 
tials eVi and eVi-, we linearize the distribution functions 
with respect to them, and obtain a linear equation for 
the potentials, which can then be solved. 

If we proceed with the analytical approximation in the 
limit A ^> Ex,eV,kBT by neglecting Tand the energy 
dependence of T>x, we obtain the dominant term: 
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Here, I S (T) = (Aa N /e) f °° jsfl& E is the observable 
equilibrium supercurrent flowing in the system when all 
parts are at the temperature T and there are no poten- 
tial differences. Moreover, f1 k = tanh (E/(2ksTk)) is 
a linearized distribution function, and R^ — L^/ (A^cr^) 
are the normal-state resistances of the wires. Thus a dif- 
ference in equilibrium supercurrents due to the varying 
fi, contributes significantly to the thermopower. 

Similarly, we can take the effect of Tinto account and 
obtain the correction terms 
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FIG. 1: N-S thermopower at different temperatures T\ and 
T2 of the normal electrodes, for a setup with Lk and Ah as- 
sumed equal for all wires, and in the limit A 3> Et, eV, UbT. 
The analytical approximation JHJ with ^correction @ devi- 
ates at worst 15% from these results, but the difference is over 
5% only where Qns is small, i.e. fcsTi,fcgT2 > 10 Et- In all 
of the figures, Et = hD/(Ls + L4 + L5) 2 , and we plot the 
results for <f> = n/2 yielding a near-maximal supercurrent. 

where Rsns = R 3 + R4 + Rs- Here, we denoted {%,) = 
■jj- J^ k T k dx (shown in Fig. 01 compared with jg). The 
correction is necessary, as it compensates for the fast de- 
cay of V®, 2 at high temperatures fcsT > 10 Et, but it is 
not negligible even at lower temperatures (see Fig - - 

The induced Qns, 2 = V 2 /(T 2 - Ti) w Qns,i is shown 
in Fig. ^ The magnitude of Qns is of the order ~ 4*- 
at highest, but there is also a strong temperature depen- 
dence. Figure |21 shows cross sections of Fig. ^ compared 
with the approximation (jSJ). 

We see that Eqs. iJSJ and @ predict an induced N-N 
potential difference 

AV = V 2 - V t w {V 2 ° - V?) + {Vi - V?) . (10) 

In a left-right symmetric structure both terms vanish, but 
Qnn may still be finite since T>t is energy-dependent: 

r°° - , /sec^f^fr) sech 2 (^V)\ 
AV w / R^ 1 y2k B nJ \2k B T 2 ) , E 

Jo V 2M 2fcsT2 / 

x (J?i + (R 3 R 5 )/(2R 3 + R 5 ))(V? + V, 1 ) , (11) 

valid for a symmetric structure R\ = R 2 , R 3 = R4. Here 
the coefficient R = R x + (R 3 R 5 )/(2R 3 + R 5 ), where 
Rk = Mx,kRk are the spectral resistances of the wires. 
The voltage (fTTJ) can be understood to be caused by the 
proximity-effect-induced temperature dependence of con- 
ductances 1 1 1| - which creates asymmetry in resistances. 
However, in reality the asymmetry of the structure causes 
likely a more significant effect (see Fig. 0}, especially at 
Ti w T 2 , where Eq. (JTTJ> predicts ~ (T 2 - T) 2 Q. 
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FIG. 2: N-S thermopower along the cross sections of Fig. 
Solid: Qns at Ti «Ta. Dotted: Approximation J^J. The cor- 
rection © accounts for most of the difference. Dashed: Qns 
at fcflTi = 3.6 Et with varying T2. Dash-dotted: the corre- 
sponding approximation. Inset: The setup under considera- 
tion: two superconducting (S) reservoirs with phase difference 
<f> connected to two normal-metal (N) reservoirs through dif- 
fusive normal-metal wires. We assume that the lengths Lk, 
k — 1,...,5 of the wires satisfy £0 J$ Lk <C 1$,Ie, where 
£0 = \/hD /(2A) is the superconducting coherence length, 
Z4, the phase-coherence length and Ie the energy-relaxation 
length. We also consider the wires as quasi- ID structures 

Equation JHJ implies that the thermopower should os- 
cillate as a function of the phase difference 4> between 
the two superconducting elements, because the equilib- 
rium supercurrent oscillates roughly as sm(^). Numerical 
simulations show (see Fig.|2t) that also the exact solution 
oscillates similarly, vanishing at = 0. 

Besides changing the prefactors in Eqs. Q8l9p . vary- 
ing the resistances of the wires from the near-symmetric 
presented in the figures changes the behavior in the su- 
percurrent and the coefficient T. In general, large 
departures from such symmetry decrease the potentials. 

A finite value for A causes two distinct modifications to 
the thermopower. First, the coefficients (JSJ are modified, 
but changes are mostly only quantitative, e.g. sharpening 
of peaks (see Fig.|3Ji). Secondly, there is also a contribu- 
tion from energies E > A, which couples the supercon- 
ductor temperature To to the system. The latter effect is 
weaker than those predicted by Eqs. (jSJ and at least 
for A > 30 Et ■ Although the coupling of To is weak, it 
induces finite potentials even for Ti = T 2 (Fig. 13b). 

Our predictions agree quantitatively with the experi- 
mental results with the correct order of magnitude for 
both the linear thermopower £| and the temperature 
scale [jjij]. The flux dependence (antisymmetric about 
<f> = — this holds for the exact as well as the ap- 
proximate solutionsjis in accord with most of the mea- 
surements in 12- 0. [j| However, we cannot explain the 
symmetric oscillations with respect to zero flux, seen in 
the "house" interferometers in Ref. 0. Moreover, the 
main result, Eq. cannot describe a sign reversal of 
Qns E| , but there is no principal reason forbidding such 
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FIG. 3: (a) N-S thermopower at fc s Ti = 3.6 E T for 
various values of A: oo (solid), 54 _Et (dotted), 27 £*r 
(dashed). Superconductor temperature To = 3.6 Et is fixed. 

(b) Potential eVi induced due to the energies above A, with 
A = 27 E T . Here, k B T = 3.6 E T but Ti = T 2 vary. 

(c) Phase oscillation of Qns (solid) and Qnn (dashed), at 
fcsTi = 4.5 Et, fcsT2 = 1.8 Et- Due to a numerical conver- 
gence problem, there is no data for \(f>\ > 2.1. (d) Spectral 
variables: Lsjs (solid), (71) (dashed), (T3) (dotted). In each 
figure, all the lengths Lk and areas Ak are assumed equal. 
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FIG. 4: Left: N-N thermopower induced by setup asymme- 
try, at fixed /cbTi = 3.6 Et- One wire (1 or 3) is assumed 
to have a different resistance than the others. Right: N-N 
thermopower with the same parameters, but with Ti ~ T2. 



an effect in a suitable structure. Nevertheless, further 
experiments are required to quantitatively demonstrate 
the connection between the thermopower and the super- 
current. 

Results of similar type as presented in this Letter for 



the N-S thermopower have been obtained for small tem- 
perature differences in Refs. 0,0, assuming high tun- 
nel barriers at the N-S contacts. However, the direct con- 
nection between the thermopower and the supercurrents 
as in Eq. (JSJ has not been shown. Moreover, Ref. dis- 
cusses a finite N-N thermopower from the energies above 
A. Our results show that for an appreciable temperature 
difference between the N reservoirs, this effect is washed 
out by the asymmetry effects, at least for A > 30Et- 

In summary, we have obtained a relation linking the 
voltages induced by a temperature difference to the 
supercurrent in a mesoscopic structure. The phase- 
oscillating N-S thermopower is mostly induced by 
the temperature dependence in the supercurrent, and 
the N-N thermopower can be attributed to left-right- 
asymmetries in a structure. These effects are indepen- 
dent of electron-hole asymmetry, and can be much larger 
in magnitude than the thermopower due to electron-hole 
symmetry breaking. 

We thank Jukka Pekola, Frank Hekking and Mikko 
Paalancn for helpful discussions. 
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